Spontaneous rotation in one-dimensional systems of cold atoms 



OO 

o 
o 

(N 



Akiyuki Tokuno 1,2 and Masahiro Sato 3 
1 Department of Physics, Tokyo Institute of Technology, 
Oh-okayama, Meguro-ku, Tokyo 152-8551, Japan 
2 Institute for Solid State Physics, University of Tokyo, Kashiwa 227-8581, Japan 
3 Condensed Matter Theory Laboratory, RIKEN, Wako, Saitama 351-0198, Japan 

(Dated: July 22, 2008) 

We theoretically study harmonically trapped one-dimensional Bose gases (e.g., 7 Li, 23 Na, 39 K, 
87 Rb, etc.) with multibands occupied, focusing on effects of higher-energy bands. Combining the 
Ginzburg-Landau theory with the bosonization techniques, we predict that the repulsive interaction 
between higher-band bosons and the quantum fluctuation can induce the ground state with a finite 
angular momentum around the trapped axis. In this state, the Z2 reflection symmetry (clockwise 
or anticlockwise rotations) is spontaneously broken. 

PACS numbers: 05.30.Jp, 03.75.Hh 
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I. INTRODUCTION 

For the last decade, cold-atom systems [J [2, B 0] have 
been vigorously studied and experimental techniques of 
controlling them have been greatly developed. As a re- 
sult, we can now tune various parameters of cold-atom 
systems, even including dimensionality. These devel- 
opments have enabled us to encounter several phenom- 
ena: For instance, superfluid-insulator transitions on op- 
tical lattices, dynamics of Bose-Einstein condensations 
(BECs), BCS-BEC crossover in Fermi gases, etc. They 
are all difficult to be realized within traditional experi- 
ments using electronic systems in solids. 

Among various new fields of cold-atom physics, in this 
paper, we focus on one-dimensional (ID) trapped Bose 
gases [E S S S 0j : For example, polarized atoms 
7 Li, 23 Na, 39 K, 87 Rb, etc. One-dimensional quantum 
systems have been shown to exhibit unusual interesting 
features due mainly to effects of strong quantum fluc- 
tuations. For example, it is widely known that BECs 
never occur even at zero temperature in ID repulsively 
interacting Bose gases in the thermodynamic limit. In- 
stead, a Tomonaga-Luttinger liquid (TLL) state is pre- 
dicted to appear. Actually TLL and Tonks-Girardeau 
(TG) gas |ilQj — the hard-core interacting limit of TLL 
— behaviors have been observed in recent experiments of 
ID cold Bose gases. [U 01] 

One-dimensional confined systems in cold-atom exper- 
iments often take a cigar shape as shown in Fig. [1] This 
shape is of course caused by a 2D (y-z plane) potential, 
which strongly restricts the in-plane dynamics and makes 
the radial energy spectrum discrete, i.e., yields a multi- 
band structure. The dimensionality of the cigar-shaped 
systems is known to be determined by the 3D s-wave 
scattering length a between atoms and the particle den- 
sity p per ID unit length. For ap <C 1 the system may 
be regarded as a purely ID gas, namely, almost all of the 
atoms are in the lowest band. On the other hand, for 
the opposite case ap > 1, 3D nature strongly survives. 
In the former case, the low-energy physics is character- 
ized by a single dimensionless parameter y -1 = aiup, 



where am is the effective ID scattering length. 0] A 
small 7 corresponds to the weakly interacting regime in 
which low-energy properties can be captured by mean- 
field (Gross-Pitaevskii-type) theory, whereas for a large 
7, effects of the interaction is considerably strong and 
thus the TLL or TG-gas behaviors clearly emerge. 

Thanks to the recent studies of purely ID atomic gases, 
their understanding has been greatly deepened. However, 
there still exist important issues in ID confined systems: 
One of them is what happens if atoms occupy not only 
lowest band but also higher ones. As yet, effects of higher 
bands have not been well discussed well. Since multi 
bands may be regarded as an internal degrees of freedom 
such as electron spins, the system is expected to contain 
rich physical properties. In this paper, we will provide 
an answer to this problem. To this end, we consider a 
simple 3D single-component Bose gas system with a re- 
pulsive interaction and a 2D harmonic potential. As we 
will explain, due to the potential, each energy band takes 
an angular momentum I: the lowest band has I — 0, the 
second lowest ones have I = ±1 and are hence doubly de- 
generate, and so on. A positive (negative) I corresponds 
to atoms rotating clockwise (anticlockwise) around the 
confinement axis. For simplicity, we will concentrate on 
the case where only the lowest band and the degenerate 




FIG. 1: One-dimensional trapped Bose gas confined along the 
x axis. In the text, we predict that the Bose gas spontaneously 
circulates along the arrow (or along the opposite direction of 
the arrow) under a certain condition. 



second lowest ones (totally three bands) are occupied. In 
this situation, it is shown that at least if the 2D har- 
monic potential is sufficiently stronger than the interac- 
tion between atoms, the particle densities of the I = ±1 
bands become imbalanced and the ground state sponta- 
neously rotates with a finite angular momentum as shown 
in Fig.[TJ We will explain in detail the mechanism of this 
spontaneous rotation phenomenon. 

This paper is organized as follows. First, we define 
the 3D Hamiltonian of the Bose atoms confined by a 2D 
harmonic potential in Sec. [TTJ Then, we reduce the 3D 
Hamiltonian to the ID effective one in order to describe 
the low-energy physics. In Sees. IHII and llVi based on the 
reduced Hamiltonian, we analyze the low-energy physics 
of the trapped bosons. Employing the Ginzburg-Landau 
(GL) theory and the bosonization techniques, we show 
that the spontaneous rotation as in Fig. [T] can occur un- 
der certain, realistic conditions. The final Section [V] is 
devoted to the summary and the discussions for our re- 
sults. In Appendix we summarize the symmetries of 
our Bose gas system, which play important roles in the 
analyses of Sees. HIT] and [TV] 



II. 



BOSE GAS AND EFFECTIVE 
HAMILTONIAN 



First, we introduce a 3D Bose gas in the presence of a 
2D harmonic potential, whose Hamiltonian is defined as 



H 



dr 



h 2 v 2 

2m 



H + V{r) ) ip + ^-V 2 + 



where ^(r) and T> = ip^-ip are the annihilation field and 
the density operator of the bosons, respectively. Positive 
parameters m, fj, and U respectively mean the mass of 
bosons, the chemical potential and the repulsive coupling 
constant. The 3D s-wave scattering length a is related to 
U via U = 4irh 2 a/m. Three or more-body interactions 
are assumed to be quite small compared with U. The 
2D confinement potential V is a harmonic type: V(r) = 
imwg(y 2 + z 2 ) — huiQ, where the axial direction is equal 
to the x axis (see Fig. [TJ . Here we note that the Hamil- 
tonian JTJ) is invariant under the reflection and the 
SO(2) rotation in the y—z plane [ip(x, y, z) — > ip(x, —y, z) 
and ip(r) — > tp( x : ycos& — z sin $, ?/ sin $ + z cos $), where 
$ € R] . Other symmetries of the model ([TJ) are discussed 
in Appendix [Al 

In order to make the one dimensionality of the 
model (JT|) more visible, it is useful to expand ip(r) 
in terms of eigenstates of the 2D harmonic oscillator 
Un,i{y, z) as follows: 



(2) 



n=Q I 



where I is the angular momentum in the y — z plane and 
can take — n, — n + 2, . . . , n— 2, n for a given n. The angu- 
lar momentum is a good quantum number because of the 




FIG. 2: Single-particle energy band structure and density dis- 
tributions (sketch) in each band. We assume that the bosons 
live only in lower three energy bands. The density p a is de- 
fined in Eq. @. 



SO (2) rotational symmetry of Eq. ([T]). The boson field 
describes the physics of the x direction. Substituting 
Eq. @ into Eq. |T]), we obtain the single-particle energy 

dispersion e n j(k) = ^^- + fku n which is shown in Fig. [2] 
It is noteworthy that all of the bands e n j with the same 
index n are degenerate, i.e., the (n + l)-fold degeneracy 
exists. For example, as we already mentioned, the second 
lowest bands ei,±i ar e doubly degenerate. 

To discuss effects of multi bands, we assume that the 
bosons occupy the lower three bands eo.o, e i,±i as m 
Fig. [2 To satisfy this situation, we impose the following 
conditions: 



Hluq < /i < 2Hujq, 
a/a± <C f (i.e. 
2Huj - n > U/a]_, 



hue > U/a\), 



(3) 
(4) 
(5) 



where the confinement radius aj_ = ^h/mujQ represents 
the width of the wave function u ny i{y, z) around the ori- 
gin (y, z) — (0,0). The first condition © ensures that 
the lower three bands have finite boson densities. The 
other ones © and mean that the energy scale be- 
tween neighboring bands is sufficiently larger than the 
strength of the interaction. These would allow us to ne- 
glect bosons in all of the higher bands e n >2,i (i-e., bands 
£n>2,; are nearly empty of bosons). 

Under these assumptions, let us continue to study the 
model ([T]). Since the low-energy physics must be gov- 
erned by the bosons on occupied bands with (n, /) = 
(0, 0), (I, ±1 ), the low-energy effective theory can be ob- 
tained by integrating out the degrees of freedom of the 
bosons on all the vacant bands. In order to carry out this 
integration, it is convenient to introduce the Euclidean 
action corresponding to the Hamiltonian (|TJ) through the 
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path integral formalism as St, 



So = I drdx 

■J ^ — n f 



n=0 I 



Sq + S\ + Si n t, 

h 



2m 



p + nhujQ 



4>n,l 



+ g^-^ [2^0,0 + Pi,i + Pi, -i + 4 Po,o(pi,i + Pi -1) 

+ 4p 1) i/>i i _i +2(00,000,001,1^1,-1 + h.c) , (6a 

/°° r h 

n=2 i 



9j, — /i + n/io^o 



S int = ^jdrd X Y:' E'«3 



(6b) 



"3,14 
i 4 



{nj ft} 



x ^mM^M^-.h^M^ ( 6c ) 

where r is the imaginary time, and p nj / = 4>Zi4>n,l is a 
boson density of a band (n,l). The prime of the sum- 
mations in Si n t means the sum over all the possible sets 
of (rii, li) (i = 1, • • • , 4) except for all of the interaction 
terms of Sq. Equations (j6a|) and ([6b]) are, respectively, 
the part including only the bosons of the occupied bands 
and the kinetic part of massive bosons of all the empty 
bands. The remaining part S- ln t contains all the interac- 
tion terms including empty-band boson fields. The cou- 



pling constants P t 



711,712,713,714 



in Si n t are expressed as 



Pi72i, 712,713,7*4 li 1 * * / —\ 

h,h,h,h = J dydZ u n 1 ,h U ni,l3 U ns,l3 , ">n*,h- ( 7 ) 

The integration of the empty-band bosons in Stot can 
be performed by making use of the representations (|5a|) - 
(|6"cl) and the cumulant expansion. The resulting action is 
written as S 1D = S + (Si^i " §((<S? n t)i ~ (^int)f) + • • • , 
where (•••)! stands for the expectation value with re- 
spect to Si. For the action Sid, effects of higher-energy 
bands are taken into account through virtual scatter- 
ing processes. In this calculation process, the Matsub- 
ara Green's function G n (r,x) = — (T T <j) ni i(x, t)0* ;(0,0)) 
of the massive bosons with e n >2,i may be approximated 
as A~ x 9 s (r)6 s(An - \x\)Q s (A n - |r|), where X n = 
y/2irh 2 /mA„ is the thermal de Broglie wave length, 
A„ = nhujQ — p is the gap of bosons, and <d s is the Heav- 
iside's step function. 

The resulting effective Hamiltonian corresponding to 
the action Sid is given as 



Hid = dx Ho + H u 



(8a) 



Ho = T^y^dx^idx^a - P0P0 - Mi (pi + P-i), (8b) 



Hint — 



2m 
U 



■ [2r p 2 Q + Vi{ P \ + p\) + 4r 0lP0 (pi + p-i) 



+ 4r±ipip_i+2r o±1 (0S^i^_i + ft.c.) , (8c) 



where a — 0,±1 correspond to band indices (n,l) — 
(0, 0), (1, ±1), respectively. The Hamiltonian Ho is the 
free boson part and po = P > and p± = p — Tiwo > 
are the chemical potentials. The first and second terms in 
Hi n t are the two-body interactions in the same band. The 
third and fourth are those between the different bands. 
The final term is the tunneling between the lowest band 
and the doubly degenerate bands. The coupling con- 
stants T Q are modified from unity (the bare value in So) 
due to vacant-band effects. 

Calculating lower-order cumulants carefully, one can 
find that, for example, massive bosons on £2,1 varies cou- 
pling constants as T a sa 1 + J2™=i Ca,n(a/a±_) n (c a>n 
are nonuniversal constants). This result strongly sup- 
ports that the cumulant expansion is reliable and works 
out under the conditions ©-([5]) and T a = 1 is semi- 
quantitatively correct. The cumulant expansion gener- 
ates other three- or more-body interactions, but they are 
also negligible compared with the two-body ones in H - m t ■ 



III. ANALYSIS I 

In the preceding section, we have obtained the ID ef- 
fective Hamiltonian ((SJ). Here, analyzing it, let us inves- 
tigate the low-energy properties of the Bose gas under 
the conditions dHJ)- (EJ) ■ First, to evaluate the mean den- 
sity profile of each band, we introduce the following GL 
potential from H m t, 



J~ GL 



-P0P0 - PlPs + 



U 



87Tfl' 



[2Topl+T s p 2 s +T a pl 



+4r 01 p oPs + 2T ± lPQ y/ ft - ft cos(20 o - 26.)] , (9) 



where T s 

Pa e , 



j = (ri ± 2r±i)/2 and we have defined <j> a — 

p s<a = p± ± p-i and 9 Si a = (01 ± 0-i)/2. 
The quantity p s is the total density in the second low- 
est bands, while p a stands for the angular momentum 
density in the same bands [the total angular momentum 
is C = 1 Jdx lp n ,i{x)\. Therefore, a solution with 
p a 7^ 0, if possible, means the emergence of the spon- 
taneous rotation and the breakdown of the Z2 reflection 
symmetry in the trapped plane. Each term in .Fgl has 
the following physical meanings: (i) Positive /io,i induce 
finite densities po,±i 7^ and the relation po 3> pi must 
yield po ^> p±i, (ii) ro, s (> 0) tend to decrease po lS , while 
r a (< 0) promotes the growth of p a , (iii) r i(> 0) favors 
the decrease of po iS , and (iv) the tunneling r ±i(> 0) 
inversely enhances po,s and decreases p a . Recalling the 
condition we can guess that bosons on the second 
lowest bands interact only with the lowest-band bosons 
with a kinetic energy eo,o(&) ~ fi&>o, and hardly influence 
the bosons with lower energies. It is hence expected that 
roi and Toil are overestimated in J-ql- Here, we simply 
regard both the density amplitude p a and the phase 9 a as 
c numbers, although for real quantum systems, we can- 
not fix them simultaneously. This approximation would 
also overestimate the effects of To±i- 
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The GL equations dTcwjdpa 
solved as 



d^GL/dOa = are 



Po 



Pa 
00 



Aiia\ r s (j, a - (2r i - 2r 0± i)/ii 
2r r s -(2r 01 -r ±i) 2 ' 
-(2r i - r ±i)^o + 2r 0; ui 



u 

47ra^_ 



u 2T Q r s - (2T 01 - r 0±1 ) 5 



0. 



tt/2. 



(10a) 

(10b) 

(10c) 
(lOd) 



Furthermore, computing d 2 J 7 Gh/dpadpf3, we find that 
the above solution is stable and minimize ^"gl- Be- 
cause of p a = 0, the GL argument shows no symme- 
try breaking. If we straightforwardly adopt the approx- 
imation T a = 1, p s becomes negative under the condi- 
tion ©. This is unphysical and attributed to large r i 
and Toil, as expected. To recover a physical solution 
Po > Ps ^ 0, we should replace r i and r ±i with small 
values, Vol and To±%, respectively: For example, when 
we set Vox = r ±i = a/a± <C 1, the GL equations offer 

a reasonable solution p s w 3U X (p,i — 5^7 A*o)- 

Let us next take into account the quantum fluctua- 
tion around the physical GL solution. To this end, the 
bosonization 0, fl3l ITij is powerful and unbiased. It 
makes the density and field operators transform as 



Pa 



a i2n(ip a —7Tp a x) 



(lla) 



n— — oo 



(lib) 



where the phase fields 9 a and ip a > {a, a' — 
0, ±1) represent the quantum fluctuation and obey 
[0 a (x),d x xp a '(x')] = [cp a (x),d x >6 a >(x')] = iir8 a ^5(x - 
x'). Here we define (p Sia = (fi ± </?-i- Substituting 
the formula (jTTJ) into Eq. ([5]), we obtain the phase-field 
Hamiltonian, which is invariant under the symmetry op- 
erations (|A1[) and (|A2[) . Remarkably, due to r a < 0, the 
coefficient of (d x <p a ) 2 becomes negative in the Hamilto- 
nian. This means that the (tp a , 9 a ) sector has an instabil- 
ity against the fluctuation of p a . It might be restored by 
sufficiently large higher-order differential terms such as 
(<9 x <£a) 4 [3) but the emergence of those terms would not 
be expected under the condition We thus conclude 
that the quantum fluctuation violates the GL solution 
and induces a ground state with a finite angular momen- 
tum p a =/= 0. It is known that similar scenarios of symme- 
try breakings can occurs in a few systems. [H M, EE El 



IV. ANALYSIS II 

In order to examine and to enhance the validity of the 
above prediction of p a ^ 0, let us reconsider the effective 
theory ([5]) using another approximation. As we already 



mentioned, Toi and Tq±\ terms are expected to be over- 
estimated in Eq. Therefore, at first we daringly drop 
the tunneling term with Tq±i and replace Toi with a 
small effective value r i. For this simplified case, the GL 
potential is written as 



■TGL = -P0P0 - PlPs + o o 

+T a pl + ATaipoPs] ■ 
The GL equations dTah/dpa = 

2na 2 , r s ^o-2foi^i - _ iira 2 ± r /Ji- 



[2T p 2 + T s p 2 s 



(12) 



lead to p = 
and p a = 0. 



u r r 3 -2r i ' Hs u r r a -2r L . 
Namely, the GL solution again suggests no angular mo- 
mentum. However, considering the stability of the so- 
lution, we see that the density profile (po, p s , p a ) corre- 
sponds to the saddle point for Tql and it is minimized 
at Pa — > io° which is unphysical. The reason for this in- 
stability is that the interband interaction T±i is stronger 
than intraband one Ti. To recover a physically proper so- 
lution, we can add the following phenomcnological term 

8^T^t-Pa ( r « 4) > 0) to f Gh . [HI It is interpreted that 
this term originates from the higher-order corrections of 
the cumulant expansion or many-body interactions ne- 
glected in the original Hamiltonian (Q]). This modifica- 
tion offers a stable, physical solution 



Pa=± tw 



(13) 



From the consideration above, we see that the GL ap- 
proach with neglecting the tunneling term naturally leads 
to a finite angular momentum. 

As in the previous analysis in Sec. IIII1 let us bosonizc 
the Hamiltonian ijHJ) based on the GL solution (|13[) . At 
this stage, we recover the tunneling term. The resulting 
phase-field Hamiltonian is represented as 

Hid « Jdx £ {K a (d x 6 a ) 2 + K-\d x ^ a ) 2 } 

a— 0.s,a 

+ 9sad x s d x 9 a + 9osd x ipod x (p s 

+ .9 cos(26> - 2<9 S ) + • • • . (14) 

Coefficients of (d x 9 ) 2 , (d x 9 s ) 2 and (d x 8 a ) 2 are propor- 
tional to po, Ps and p s , respectively, while g sa oc p a . 

Therefore, when p a is smaller enough than po,s> any 
instability of restoring p a — does not originate from 
the differential terms. In addition, the symmetry opera- 
tions (|A2p indicate that any vertex term with ip a or 9 a is 
forbidden in the Hamiltonian (fl"4|) . (l8| Thus, we can say 
that p a 7^ is not destroyed by the quantum fluctuation 
and it agrees with our previous approach in Sec. lIIII The 
symmetries also tell us that for all of the vertex operators, 
only cos[2n((9o — 9 S )] (n € Z) are allowed to appear in 
the phase-field Hamiltonian. The final term of Eq. (|14p is 
the most relevant in the allowed vertex interactions. Gen- 
erally, for ID boson systems with repulsive short-range 
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TLL with a finite 
angular momentum 
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FIG. 3: Predicted ground-state phase diagram of the 
model (UJ) under the condition ((4]). Quantum phase transi- 
tions are expected to occur around the black points. They, 
however, might be interrupted as a crossover phenomenon 
rather than a transition, since boson densities of higher bands 
would be nonzero even when the chemical potential /i is 
smaller than huo- When p is larger than nhuio (n > 2), the 
bosons on higher bands e n >2,i would also possess an angu- 
lar momentum. It is therefore expected that, for example, 
the bosons on £2,! rotate clockwise, while those on £3^ do 
anticlockwise. Namely, the almost complete or partial cancel- 
lation of the angular momentum can occur. 



interactions, the TLL parameter -Ko,s,a are known to be 
larger than unity. [9( The vertex term cos(2#o — 29 s ) thus 
must be strongly relevant [HI and 9q — 9 S is locked at 
7r/2, which is consistent with our previous result (|10d[) . 
Introducing new fields {ip± , 6± ) = (lpq±(p Si 9q±9 s )/ y/2, 
we can interpret that the cosine term opens a gap in the 
(tp-,9-) sector. Following the method in Ref. [20|, we 
can trace out the degrees of freedom of the gapped sec- 
tor. After that, we finally obtain a two-component TLL 
Hamiltonian with the (<p+, 6+) and (ip a , 9 a ) sectors as the 
low-energy effective theory. 

On the other hand, when p a is too large in Eq. (flU|) , we 
should adopt (pi,p_i) = (p s ,0) or (0,p s ) as the proper 
GL solution instead of Eq. (fT5|) . For this case, neglecting 
all of the parts with the field </>_i in the Hamiltonian ||5J), 
we again derive a two-component TLL Hamiltonian. 

Consequently, we conclude from the second analysis 
above that the ground state possesses a finite angular 
momentum p a 7^ and the low- energy excitations are 
described by a two-component TLL. [211 ] 



V. SUMMARY AND DISCUSSIONS 

In this paper, we have studied the ID harmonically 
trapped Bose gas ([T]) with bosons on multi-transverse 
modes. Starting from the 3D boson systems with 2D 
confinement potential, the ID effective Hamiltonian is de- 
rived by integrating out the degrees of freedom of all the 
vacant bands. Applying two theoretical ways based on 
the GL approach and the bosonization to the ID effective 
Hamiltonian l[5]). we have shown that when bosons are 
filled in the lower three bands under the conditions ([3])- 
([5]), a angular momentum becomes finite and the Z2 re- 
flection symmetry in the trapped plane is spontaneously 
broken. This phenomenon essentially originates from the 
larger repulsive interaction between the degenerate bands 
(n, I) = (1, ±1) than the intraband repulsion. The second 
analysis in Sec. IIVI predicts that the low-energy excita- 
tions on the rotating ground state are described by a two- 



component TLL. [2l| It is known that when the bosons 
occupy only the lowest bands, the low-energy physics is 
governed by a one-component TLL. We thus can draw 
the phase diagram in Fig. [3] 

The conjugate field for p a is not realistic in the 
model (HJ). Any trigger of the rotating ground state 
therefore seems to be absent. However, the real trap 
potential must deviate from the harmonic type. Such 
a small deviation would help the system circulate clock- 
wise or anticlockwise. Thermal fluctuations and large de- 
formations of the trapped potential would generate kink 
structures between clockwise and anticlockwise rotating 
regimes. When Bose atoms carry an electric charge, the 
predicted rotation means spontaneous loop current and 
magnetic flux. If we suddenly change the form of the trap 
potential and then observe the real-space boson density 
profile, we can detect, in principle, whether the bosons 
rotate clockwise or anticlockwise. 



APPENDIX A: SYMMETRY ARGUMENT 

Here, we briefly summarized the symmetries of our 
Bose gas system ([T]) . They are often used in our analyses 
of Sees. HI] and [IVl 

The Hamiltonian JTJ) is invariant under the following 
symmetry operations: The translation and the parity 
for the x-axis [ip(x,y,z) — > ip(x + 8,y,z) (S € R) and 
■0(x, y, z) — ► ip(—x, y, z)], the Z2 reflection and the SO(2) 
rotation in the y-z plane [ip(x, y, z) — ► i/j(x, — y, z) and 
ip(r) — > ip(x, y cos $ — z sin <£>, y sin <E> + z cos $) ($ 6 R)], 
and the global U(l) gauge transformation [ip(x,y,z) — > 
^(x,y,z)e ie (9g R)}. 

For the bosonized phase-field Hamiltonians (e.g., 
Eq. (fT4)) ). it is important to note that all of the above 
symmetry operations can be translated to the bosoniza- 
tion language. [l8| The translation and the parity opera- 
tion for the x-direction, the Z2 reflection and the SO (2) 
rotation in the y-z plane, and the global U(l) gauge 
transformation are, respectively, expressed as 



[<Pa(x), 
[<Pa(x), 

[<Pa(x), 
[<Pa(x), 
[fa(x), 

where a 
form Eq. < 

[<Pa(x),< 
[<Pa{x),t 
[<Pa{x),\ 
[<Pa(x),t 



9 a (x)] —> [tp a (x + 6) + np a S,9 a (x + 6)], (Ala) 

9 a {x)] -> [-<p a (-x),0 a (-x)], (Alb) 

9 a {x),p a ] -> [ip- a (x),9- a (x), p^ a ], (Ale) 

)} -> [(p a (x),9 a (x) - (Aid) 

)] — * [v>a{x),6 a (x) + 9], (Ale) 

= 0, ±1. Using another set of the phase fields 



Alii as 



— ), we can further trans- 



* a (x)] -> [p a (x + 5) + -Kp a 8,0 a (x + 5% (A2a) 

t a (x)]->[-tp a (-x),0 a (-x)], (A2b) 

'a{x),pa] -» [-(p a (x) 7 -9 a (x),-p a ], (A2c) 

) a (x)]^[<p a {x),9 a {x)-$], (A2d) 

) a (x)]^[ip a (x),9 a {x)], (A2e) 
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for ((fa,0a 


), and 








r / \ 
[<Ps{x), 


9 s {x)\ - 


r/i / i 

* [0 8 {x J t 


-d),(/? s (a; + dj- 


1- 7rp s d], (A3a) 


[<Ps(x), 


9 s (x)] - 


* [-<Pa{- 


-x),0 a (-x)], 


(A3b) 


[<Ps{x), 


Ux)] ~ 


* [<p s (x) 




(A3c) 


[<Ps{x), 


9s(x)] - 


* [<Ps{x) 




(A3d) 


[<Ps[x), 


Ux)} - 






(A3e) 



for (ip s , 8 S ). Note that p a changes the sign under the Z-i 
reflection in the y-z plane shown in Eq. (|A2c[) . 

For instance, applying these symmetries, we can 



strongly restrict possible operators in the bosonized 
Hamiltonian (fT4|) . 
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